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^ \ Abstract 

rH \ 

We investigate the role the interweaving of surface vibrations and nu- 
. ^ cleon motion has on Cooper pair formation in spherical superfluid nuclei. 

^ ■ A quantitative calculation of the state-dependent pairing gap requires to go 

beyond the quasiparticle approximation, treating in detail the breaking of 
the single-particle strength and of the associated poles. This is done solving 
self-consistently the Dyson equation, including both a bare nucleon-nucleon 
interaction (which for simplicity we choose as a monopole-pairing force of con- 
stant matrix elements g) and an induced interaction arizing from the exchange 
of vibrations (calculated microscopically in QRPA) between pairs of nucleons 
moving in time reversal states. Both the normal and anomalous density Green 
functions are included, thus treating self-energy and pairing processes on equal 
footing. We apply the formalism to the superfluid nucleus 120 Sn. Adjusting 
the value of g so as to reproduce, for levels close to the Fermi level, the em- 
pirical odd-even mass difference (A ~ 1.4 MeV), it is found that the pairing 
gap receives about equal contributions from the monopole-pairing force and 
from the induced interaction. This result is also reflected in the fact that if 
one were to reproduce the observed A allowing the nucleons to interact only 
through the bare monopole-pairing force, a value of g ~ 0.233 MeV (~ 28/A 
MeV) is needed, 50% larger than the value g « 0.166 MeV (« 20 /A MeV) 
needed in the full calculation. To keep in mind that the bare and the induced 
pairing contributions to A enter the corresponding equations in a very non- 
linear fashion. It is furthermore found that self-energy processes reduce the 
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contribution of the phonon induced interaction to the pairing gap by ~ 20% 
as compared to the value obtained by only phonon exchange without taking 
into account the breaking of the single-particle strength. 

PACS codes: 21.30.Fe, 27.60.+j 

Keywords: Dyson equation, pairing gap, phonon-induced interaction 

1 Introduction 

It is well known that single-particle motion in atomic nuclei is strongly renormal- 
ized by the coupling to low-lying surface vibrations, which leads to the fragmentation 
of the single-particle strength and affects basic quantities like the level density at 
the Fermi surface, the effective mass and the width of giant resonances [1-3]. Re- 
cently, it has been shown that the induced interaction resulting from such a coupling 
(phonon-induced interaction) leads to pairing gaps which account for about half of 
the experimental values [4] . In what follows we shall refer to these calculations of the 
state-dependent pairing gap as the Bardeen-Cooper-Schrieffer plus Bloch-Horowitz 
(BCS+BH) approximation, in keeping with the fact that this form of perturbation 
theory was used to calculate the induced interaction entering the gap equation and 
arising from the exchange of phonons between pairs of nucleons moving in time re- 
versal states close to the Fermi level. In these calculations, empirical single-particle 
levels were used, and the interplay between core polarization and self-energy contri- 
butions, as well as the energy-dependence of the effective mass [2] was neglected. 

In the present paper we take a step towards a more consistent calculation of the 
renormalization effects associated with the particle-vibration coupling phenomenon, 
treating self-energy and induced pairing interaction processes on equal footing (cf. 
Fig. 1(a)), by solving the Dyson equation (also called Nambu-Gorkov equation, cf. 
e.g. [5]), describing the motion of nucleons and their contribution to both normal 
and anomalous densities. Moreover and following ref. [|j], we add to this interaction 
a monopole-pairing force acting on the orbitals close to the Fermi level. 

No vertex corrections (cf. e.g. graph (b) of Fig. 1 ) have been considered in 
the present calculations, because in the case of 120 Sn, the nucleus we use in the 
present paper to exemplify the theory, vertex correction contributions are found to 
change (reduce) the value of the state-dependent pairing gap around Fermi level 
only modestly.0 

The Dyson equation approach has been customarily used in condensed matter 
physics to deal with strong coupling superconductors (cf. e.g. ref. [5-9] and refs. 
therein). In nuclear physics, this scheme has also been used previously both in the 

1 One could argue concerning the validity of the approximation to neglect the vertex correction 
in terms of Migdal theorem [6,7] which states that the contribution of the vertex renormalization 
graphs to the total electron-phonon coupling strength is smaller than the lowest-order strength by 
a factor MI/Ef, where TiQ is a typical phonon energy, Ef being the Fermi energy measured from 
the bottom of the single-particle potential. Because in the nuclear case Ep ~ 36 MeV and frfl is of 
the order of few MeV, the ratio above turns out to be approximately equal to 10 _1 . However, this 
estimate can hardly be considered quantitative, because of the marked shell structure displayed by 
the single-particle levels, and by the associated particle-vibration coupling matrix elements (spin- 
flip vs. non spin-flip processes) as well as the strong cancellation existing between effective mass 
and vertex correction contributions associated with core polarization phenomena (cf. e.g. [3]). 



2 




(i n' of 



b) 




<- 



<- 



Figure 1. a) The self-energy ?iSP ho (w) used in the calculations. The 
wavy line denotes the unperturbed phonon Green function. The heavy 
arrowed line indicates the perturbed nucleon Green function. The small 
open circle is the vertex of the particle-phonon interaction, b) Example 
of self-energy diagram which includes the vertex correction. 



case of nuclear matter and of finite nuclei [10-15], the work of ref. [15] being close 
in spirit to the one discussed here. However in this reference the focus was set 
on applying the method to correct the phenomenological values of single-particle 
energies and of the pairing gap to avoid double counting arising in connection with 
particle-vibration coupling processes. Furthermore, since we do not make use of 
the linear approximation adopted in that reference to solve the Dyson equation, the 
present approach should be also able to deal with strong-coupling situations. 

In section 2 we present the formalism used in solving the Dyson equation. In 
section 3 we discuss the results obtained including only the phonon-induced inter- 
action, while in section 4 we add the monopole-pairing force to it. Conclusions are 
collected in section 5. 

2 Solution of the Dyson equation 

The core of the present study is the Dyson equation, which for paired systems is 
written in a matrix form, according to Nambu-Gorkov theory [5]: 



uj being the energy variable. G^iyj) denotes the perturbed single-particle (quasi- 
particle) Green function, and its major poles yield the quasiparticle energies of the 
system. The diagonal and off-diagonal elements of Eq. (2) are the particle-hole and 
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the pairing Green function respectively. G^Juj) is the unperturbed single-particle 
Green function, with elements given by, 

tG! U M = r — , (3) 

fl UJ - (£° - S F ) + ITJf, 



-< 2 H = -g° 21 m = o, (4) 

1 oo p-iV" 

n * uj + (e» - e F ) - vq^ 



where e° denotes the unperturbed single-particle energy of the single-particle state 
with quantum numbers fi ( = nlj). The Fermi level measured from the top of the 
potential is indicated by Bp. The parameter rj^ is defined as 

[ V > £ fj, > £ F > 

where rj is real and positive. The self-energy term of the phonon-induced interaction 
(cf. Fig. 1(a)) is defined as B, 



S llpho M E 12p h o (a;) 

00 du' 1 x ^ h 



X 



-oo 27r ^ h Xn 2n Xn B\ n 2j M + 1 

! ' DKu-u'), (7) 



where r 3 = ( J -l ) ■ ener gy an d the inertial parameter of the phonons are 
denoted by hti\ n and by i?An, respectively, with the multipolarity A and an additional 
index n. The particle-vibration vertex is given by 



11 -(/ipo^nii//) , (8) 



2tt Xn B Xn dr 



where i?o an d U(r) are the nuclear radius and potential. The phonons are not dressed 
and are calculated within the QRPA, as was done in ref. The corresponding 
(unperturbed) phonon propagator is 



h " UJ — UJ 1 — flQxn + —UJ + Ci/ — ftf^An + 

where 770 stands for a real positive parameter. 

We solve Eqs. (1) and (7) self-consistently, starting by inserting at the place of 
G^i in Eq. (7) the BCS Green function in an analytical form (see Chap. 7 in ||), 
calculating the integral of the self-energy TP^°{uj) analytically. The poles ±^g+ ( we 
set ReujQ^ > 0, m being another label distinguishing the pairs of the poles ) of the 
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new Green function are determined by searching, on the real axis, for the roots of 
the equation^ 

detG~\±Reu>%™) = , (9) 

where G (u) is calculated from Eq. (0) but neglecting the non-hermitian compo- 
nents: 



UJ 



-hE 12pho (uj) 



(10) 



UJ 



with e ^ = e ^ — e f . A very small value (1 keV) of the averaging parameters r\ and t)d 
is used with an idea to deal with stationary states. The imaginary parameters are 



so small that we can identify the residue of G^(u)/Ti with that of G l }-{u))/%, that 



is. 



R}L(±uj, 



"lim 



G+) 



d h 



'ID 



±Re <^q™ 



where the derivative is taken also on the real axis. i?*^(±c<;G.™) is obtained through 
a similar equation. The imaginary part of the poles ±Wg™ is fixed to be =F 1 keV 
for simplicity. Eqs. (|9|)-(|ll"|) are generalizations of those found in Sections 3.4 and 
4.3.6 of ref. Q. The new Green function can now be written as 



' Dll (. 



UJ — UJ, 
>12 



fj,m 
G+ 
fim\ 



+ 



Dll / , 

n fim[~ UJ G+J 
, am 
UJ + UJq + 



>12 



R^m( UJ G+) 

n.m. "I 

UJ — UJ', 



G+J 



G+ 



UJ + UJ, 



/im 
G+ 



(12) 
(13) 



Assuming that the ground state of the system is time reversal invariant, the quanti- 
ties G 22 {uj) and G 2 Hu) can be written in a similar way in terms of R^(±Uq!{_) and 
Rj? m * {jLuJq^) , respectively. A new self-energy is then generated, using Eqs. (12) and 
(13), as 



TiE 



llpho. 



UJ) 



plm An 4 An- An 



^II^^All/i') 



1 



X 



Dll f , ,A»'JTl\ 



2ju + 1 



UJ + UJ&™ + ftfi A n 



ir]D 



UJ — UJ 



G+ 



TiVLxn + ir]D 



.(14) 



The other elements can be obtained in a similar way, and the process is iterated until 
convergence. The Fermi level Bp is fixed so as to obtain the desired expectation value 
of the nucleon number 



(15) 



fj,m 



As can be seen From Fig. 1(a) and Eq. (|I]), this iterative process involves the 
coupling to an increasingly larger number of phonons at each iteration step. This 

2 We use a mesh of 1 keV, and once the existence of a root is detected in an interval, the bisection 
method is applied. 
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in turn enhances the fragmentation of the single-particle strength, leading to an 
increasingly larger number of poles of the Green function until saturation.^ For each 
of the poles ±Uq^, the sum R^{wg+) + R]} m (~ u G+) ls interpreted as the single- 
particle strength of the pole ( RV m {—^>G+) = ^m(^G+) ). As a rule, and depending 
on whether pairing correlations are important or not, one or two poles close to Sp 
carry most of the single-particle strength in the case of orbit als close to the Fermi 
level. In particular, in the BCS approximation the two poles ±c<jg+ carry the whole 
strength, and Rf(u& + ) + flj 1 (-£<;£+) = wj + < = 1, and \R*?(±u&+)\ = u^. In 
general, however, one needs to take into account more poles, in order to exhaust a 
substantial fraction of the sum rule 

EKmK+)+4-(-^+))=l- (16) 
m 

We shall label with //mo the pole carrying the major single-particle strength, and 
call it the quasiparticle pole. 



3 Induced interaction 

We have applied the formalism described above to the calculation of the neutron 
Green functions of the nucleus 120 Sn. In this section we shall concentrate on the 
phonon-induced interaction, while in the next we shall add the monopole-pairing 
force with constant matrix elements as a model of the bare neucleon-nucleon force. 

The unperturbed single-particle basis e° has been calculated with a Woods- 
Saxon potential using an effective mass rrik = 0.7m, m being the bare nucleon mass. 
With this value, self-energy effects are expected to lead to a sensible single-particle 
density close to the Fermi level (cf. also Section 4.6.3 in 0]). We have included in 
the calculation all the single-particle bound levels. 

The computation time needed for our calculations depends strongly on the num- 
ber of phonon modes An included. The full QRPA response for the multipolarities 
A 71 " = 2 + , 3 _ , 4 + , 5~ in the energy interval 0-20 MeV used in ref. [4] consists of about 
two hundreds phonon modes of energies hfl\ n and zero-point amplitudes /3\ n for each 
multipolarity. We include the four lowest phonons, one for each multipolarity, which 
give the largest contributions to the induced phonon interaction. We account for 
the effects of the other roots including only a few effective phonons of energy hQ c ^, 
distributed in the interval 0-20 MeV, choosing their effective strength so that when 
they are used in the calculation of ref. [4] they reproduce the state-dependent gap 
obtained there. This is obtained, considering that the sum of the (asymmetrized) 
matrix elements of the induced interaction between two pairs (jz,)j =0 , Cv)j =0 due 
to the phonons lying in an energy interval [Q a ,f2&], calculated according to the BH 
formalism, is given by 

3 The number of poles saturates because r\ and rjo are finite, and we use the mesh method for 
finding the poles as explained before. That is, the poles with very small residue (strength) are 
neglected. 
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Table 1. The energies of the phonon modes hO,^ and coupling strength 
P&/V2X + T. Tables a), b), c) and d) are for A 71 " = 2 + , 3~, 4+ and 5", 
respectively. The lowest-energy modes (n = 1) were taken from a QRPA 
calculation directly. The coupling strengths as well as the energies of the 
other modes were determined by the procedure shown in Eqs. (|l7|)-(18). 
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0.0365 
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0.0300 
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13.0 


0.0166 
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15.0 


0.0270 
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21.0 


0.0232 



where E and e v are the correlation energy of the ground state and the absolute 
value of the single-particle energy with respect to the Fermi level, respectively. The 
effective strength of the phonon representing this interval is then chosen so as to 
satisfy the equations 

(f)eE\2 n b n 2 



E - (e„ + e v < + Mlf n ) n E - (e„ + e v > + htt Xn ) 



-veff 



nnf n = nn b , (i8b) 

for the pairs (j u )j=o, Cv)j=o giving the largest contribution to the pairing gap for 
the multipolarity A. The energies and zero-point amplitudes devided by \/2A + 1 of 
the effective phonons are listed in Table 1. Then the coupling strength 



2ngB& 2A + 1 



is used for Eq. (|l4|) . The BCS+BH calculation performed with this restricted en- 
semble of phonons reproduces the state-dependent pairing gaps of ref. [4] within a 
few per cent. 
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Figure 2. a) The spectral function of 2d-^/2 calculated using the poles 
and residues of the solution of the Dyson equation. To display the results 
in terms of a smooth continuous curve, the results have been folded with 
a smooth function with FWHM of 0.1 MeV. Similar foldings have been 
performed in connection with the other plots, b) Same but for the 2(i 5 / 2 
orbital , c) The spectral function of the l/tii/2 orbital, d) The even self- 
energy of the l/tn/2 orbital (real part), e) The pairing spectral function 
calculated in a similar way as the spectral function, f ) The pairing self- 
energy (real part). 



Because the number of the poles can vary from one iteration step to the next, 
the solution displays (small) fluctuations. We obtained a reasonable accuracy in the 
particle number (N) = 70 ± 0.1 and negligible fluctuations in the perturbed single- 
particle levels and the state-dependent pairing gaps (definitions are given later) for 
most of the orbitals. The sum of the single-particle strength for a given orbital /x is 
rule larger than 0.90 in the valence shellfl 

4 To be noted that of all the (hundreds or eventually thousands) poles, only a fraction (< 
10—20) carry an appreciable single-particle strength for orbitals near the Fermi level. Consequently, 
calculations made considering only these poles explicitly lead to results which coincide within 10%, 
with the results of the more accurate calculation presented in this paper. 
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The spectral function 



S,(u) = \\lmGf(uj)\ , (19) 
irn 



of some levels is shown in Fig. 2 (graphs (a),(b) and (c)). From this figure it can be 
seen that the quasiparticle approximation (keeping only one (no pairing) or two (with 
pairing) pole(s)) works well for the level 2c? 3 / 2 , (particle-like quasiparticle, cf. Fig. 
2(a)), and lh n / 2 (with an associated single-particle strength R l } mo {~ Uq^ ) ~ 0.3, 
cf. Fig. 2(c)), but breaks down for the hole-like 2<i 5 / 2 level (cf. Fig. 2(b)), which 
displays a substantial degree of fragmentation, because one of the poles of 2<i 5 / 2 
is almost degenerate with the lhu/ 2 X 3~ state, then the single-particle strength 
(residue) of the pole becomes small (see Eqs. (0), (0) and ([□])), and the other 
poles carry more strength. In addition, a large matrix element connecting lhn/ 2 
and 2c?5/2 with X w n = 3~ 1 (see Eq. (|8|)) enhances the fragmentation. Besides the 
spectral function, we also display in Fig. 2 the even self-energy 

^llphoeven^ = |^llpho (w) + E^-^)) ■ (20) 

This is shown in Fig. 2(d) for /i = lhu/2- Its value calculated at the quasiparticle 
peak of the spectral function gives the approximate energy of the renormalized 
single-particle spectra 

4 = ttt^ + hRe £^ h ° eve >£7)) + e F , (21) 



with 



KRz £!, lph ' 



D odd/ fim 



The function hY,j} pho ° dd (uJq^) is the odd component of the self-energy. (See Chap. 
7 in ref. M in connection to Z^uJq^).) 

The original single-particle spectrum is compared in Fig. 3(a) with the energies 
e 1 ^ of the main (quasiparticle) peaks. It is seen that the self-energy processes increase 
the level density in the valence shell. This increase is essentially due to the rise of 
the energy of the hole states. 

The spectral function and self-energy shown above have their counterpart in the 
pairing channel. The pairing spectral function Sjj 2 ^) corresponds to the anomalous 
density. In the BCS limit it is peaked at the quasiparticle energy and corresponds 
to the product u^v^. The function S^{ui) is shown in Fig. 2(e) for \i = lh u / 2 . The 
pairing self-energy (u;) is shown in Fig. 2(f), again for the level /i = lhn/ 2 . 

The complex state-dependent gaps are given by [5] 



ft E 12pho, 

a pho _ v o+ / (22) 

H 7 (, ,/«no\ - \ I 

Z A W G+ ) 

The assignment of a state-dependent pairing gap as well as a perturbed single- 
particle energy is valid only to the extent that the quasiparticle approximation is 
valid, thus the assignment may not be adequate for the single-particle orbital d 5 / 2 
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Figure 3. a) The unperturbed single-particle levels e° near the Fermi 
level ef and the perturbed levels corresponding to the quasiparti- 
cle poles, b) The state-dependent pairing gaps corresponding to the 
quasiparticle poles (open circles) and the gaps of the BCS+BH cal- 
culation (filled circles) as a function of the renormalized (perturbed) 
single-particle energies c) A magnification of (b) near the Fermi 
level ef- The arrow indicates the Fermi level associated with the Dyson 
calculation. 
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(cf. Fig. 2(b)), which is strongly fragmented. In this tentative quasiparticle 

pole (the right major peak in Fig. 2(b)) was used. The results of the calculation 
are collected in Table 2, and the gaps associated to the quasiparticle peaks oJq^ 
are shown in Figs. 3(b) and (c) in absolute value, where they are compared with 
the pairing gaps in the BCS+BH approximation calculated using Table 1 and the 
renormalized (perturbed) single-particle spectra e 1 ^ of the Dyson calculation (these 
energies were used for calculating e v in Eqs. ( pf ) and ( |18a| ). Close to the Fermi level, 



see Fig. 3(c), the gaps of Dyson calculation are systematically lower; differences of 
up to 400 keV between the two calculations are observed. In average in the valence 
shell, the difference is ~ 200 keV. 

It is also interesting to calculate more global quantities, like e.g. the pairing 
correlation energy. Within the framework of the Dyson equation, the contribution 
to the total ground state energy of the system arising from pairing correlations is 

<° = ~i E(2i, + 1)^ ^ (%^{u)\Gf{u)^ 

+ ^^ h »lG; 2 ( W )e" v ) , (23) 

in keeping with the fact that the quantities entering the above equation are the 
only ones depending explicitly on the anomalous density. In the case of 120 Sn, 
= -3.9 MeV, a number to be compared with £ pair (BCS + BH) = -4.6 MeV 
calculated making use of the BCS relation 

£ pair (BCS + BH) = -J2 ^-U^A^BCS + BH) , (24) 

where A M (BCS + BH) is the pairing gap calculated in the BCS+BH approximation. 

We have also performed calculations using only two phonons, corresponding to 
the low-lying X n = 2 + and 3~ modes (n — 1 in Table 1). The resulting average value 
of the pairing gap in the valence shell is 0.29 MeV, which represents about 50% of 
the value (0.54 MeV) obtained using all the effective phonons listed in Table 1. It is 
seen that, while the low-lying phonons give the largest contribution to the induced 
pairing gap, the inclusion of the other modes increases the gap further. 

The perturbed single-particle levels in the valence shell in the calculation using 
only the low-lying 2 + and 3~ modes are less compressed around the Fermi level by 
only 100 keV in average as compared to those obtained in the full calculation, thus 
the perturbed spectrum is similar to Fig. 3(a). 



4 Induced plus monopole interaction 

The results presented in the previous section included only the contribution of the 
phonon induced interaction. However, a realistic calculation of pairing correlations 
must include also a bare nucleon-nucleon interaction. We do this in the present 
section, adding to the self-energy of the phonon induced interaction used above (cf. 
Eq. (7)) the standard monopole (seniority) pairing field acting in the valence shell 
(v.s.). Associated self-energy reads 



n 
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Figure 4. a) The pairing gaps obtained solving the Dyson equation in- 
cluding the monopole+induced interaction (open circles) are compared 
with the corresponding BCS+BH calculation (filled circles), b) A mag- 
nification of (a) near the Fermi level indicated by the arrow. 




/i'ev.s.,m (2t> ) 

for other \i , 



ftE^ u = ft£^ u = . (26) 

The experimental value of the pairing gap, deduced from the experimental odd-even 
mass difference, is A ~ 1.4 MeV (cf. Chap. 2 in WW )- The state-dependent gap 
obtained using g = 20 /A = 0.166 MeV, plus the induced phonon interaction used in 
the previous section, is shown in Fig. 4; the average pairing gap in the valence shell 
1.46 MeV is close to the experimental value. The strength g = 0.166 MeV should be 
compared with the value g = 28 /A = 0.233 MeV needed to obtain the experimental 



value using the monopole-pairing field ( p5f) and (26) without any phonon coupling 
(that is, performing a simple BCS calculation starting from the unperturbed single- 
particle spectrum with m* = 0.7m). 
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The enhancement of the pairing gap due to the induced phonon interaction 
seems to be in contrast with the effect of the induced interaction in nuclear matter, 
where it usually decreases the pairing gapQ. In fact, in neutron matter the induced 
interaction takes place mostly through the exchange of 5 = 1 excitations via the 
spin-spin part, (s 1 ■ s 2 )V s , of the nucleon-nucleon force, which leads to a repulsive 
interaction, and tends to suppress the pairing gap. In the case of finite nuclei the 
absence of collectivity associated to the 5 = 1 modes should make its contribution 
only marginal with respect to the effect of the strongly collective low-lying surface 
modes treated in this paper. 

In order to get more insight in the result shown in Fig. 4, we decompose the 
anomalous self-energy £* 2 in the sum of two terms, including only the matrix ele- 
ments of either the induced or the monopole interaction, 

KEfiug?) = hEf pho {uj^l ) + hZf P0 . (27) 

The imaginary parts of fr£^ 2pho (u;g™ ) and hT}^ P0 are negligible, and the real parts 
have the same sign. Also the state-dependent pairing gap, defined as 

a l£ v (j + (9£l\ 

Z 'H\ U G+ ) 

can then be separated into two contributions, 

A, = Af ° + A™ , (29) 

with 

fcyl2P0 



The monopole interaction has only a small effect on the values of Z M , since the 
interaction does not act in the particle-hole channel: from Tables 2 and 3 it can be 
seen that the values of Z^ obtained with or without the monopole interaction agree 
within 10% for all the valence orbitals. From Table 3 it can also be seen that the 
pairing gaps A M receive similar contributions from A^ 110 and from A P0 , except for 
the fragmented level <i 5 / 2 where the phonon gap is much larger. In particular, the 
values of Ajf 10 are enhanced compared to those obtained in section 3 without the 
monopole interaction, and listed in Table 2. 

Finally, in Fig. 4 we compare the gaps obtained in the present calculations with a 
BCS+BH calculation of the kind performed in [4], but including the same monopole 
interaction and the same effective particle-vibration couplings as used here. It is 
seen that the Dyson calculation leads to gaps which are lower again by about 200 
keV in average in the valence shell as compared to the BCS+BH results. This is 
because nucleons spend, in the Dyson scheme, part of their time in more complicated 
configurations. 



5 It was shown in ref. |l6[ that the sign of the polarization effect is density dependent in nuclear 
matter and that the effect reduces the gap at low densities. 
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Table 2. The result for the orbitals in the valence shell of the Dyson 
calculation with only the phonon-induced interaction including the ef- 
fective coupling strengths. Only the real parts of the gaps are shown, 
since the imaginary parts are negligible. fiS^ 2pho denotes the value at 
the quasiparticle pole hT}^°{u^l ). 



[MeV] [MeV] [MeV] 

-9.232 -0.642 -1.049 1.63 

13.166 -0.951 -2.396 2.52 

-9.977 -0.519 -0.866 1.67 

11.618 -0.183 -0.231 1.26 

-9.000 -0.353 -0.578 1.64 



hn/2 
d 5 /2 - 

Sl/2 
S7/2 
d 3/2 



Table 3. The result of the Dyson calculation with both the phonon- 
induced interaction, including the effective coupling strengths, and the 
monopole (seniority) pairing force. For the notations, see the text. 
Note that KT, 1 ^ , and Z M are abbreviations of ^S^ 2 P ho (a;^ ), and 
Z/j,{oOq^°), respectively. Again we show the real parts of the pairing 
gaps. 



[MeV] [MeV] [MeV] [MeV] [MeV] 

1.400 -0.868 -1.287 -0.531 -0.788 1.48 

1.942 -1.602 -3.710 -0.340 -0.788 2.32 

1.429 -0.919 -1.419 -0.510 -0.788 1.54 

1.244 -0.677 -0.942 -0.566 -0.788 1.39 

1.311 -0.792 -1.204 -0.518 -0.788 1.52 



hn/2 - 
d 5 /2 

s l/2 
§7/2 

d 3/2 - 



5 Conclusions 

One can conclude, that in the case studied the induced phonon interaction con- 
tributes substantially to pairing correlations: the strength of the monopole interac- 
tion, needed to obtain the experimental value of the pairing gap, is strongly reduced 
when the phonon contribution is included, as the gaps receive almost equal contri- 
butions from the monopole and from the induced interaction. 

It is also seen that a quantitative estimate of the pairing induced interaction re- 
quires a self-consistent treatment of the particle- vibration renormalization processes, 
as the one performed here solving the Dyson equation. In fact the self-consistent 
treatment lowers the gaps by about 20% as compared to the empirical treatment 
of the effective mass processes, which does not include a detailed description of the 
breaking of single-particle strength. 
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